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r£^ • Abstract 

■*— > ■ 

We derive a collection of identities for bivariate Fibonacci and Lu- 
cas polynomials using essentially a matrix approach as well as prop- 
erties of such polynomials when the variables x and y are replaced by 
polynomials. A wealth of combinatorial identities can be obtained for 
£> ' selected values of the variables. 
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We define bivariate Fibonacci polynomials as 

Fn(x, y) = xF n _i(x, y) + yF n _ 2 (x, y), F (x, y) = 0, F t (x, y) = 1, 



and bivariate Lucas polynomials as 



L n (x, y) = xL n _i(x, y) + yL n ^ 2 (x, y), L (x, y) = 2, Li(x, y) = x. 



e 
> 

S^ . We assume x ^ 0, y ^ 0, x + Ay f: 0. 

The roots of the characteristic equation are 



X + y/x 2 + Ay X - yjx 2 + Ay 

a = a{x,y) = , = /3(x, y) = . 



We have a + j3 = x, a/3 = —y and a — j3 = y 'x 2 + Ay. The Binet's forms are 

F n (x, y) = ?—^, L n (x, y) = a n + (3 n . 
a — p 

Many basic facts concerning these kinds of polynomials can be found in [1]. 



2 Some Matrix Techniques 

Let 



B 



Then 
where 

It follows 



x + 2y x 

xy 2y 

B = 2yl + xA, 



x 1 

y o 



B" = ^ (2y) n - k x k A h 

It is known that (e.g. [1]) 



F n+1 (x, y) F n (x, y) 
yF n (x, y) yF n _i(x, y) 

Now using the method exposed in [2] to B we find, after some manipulations, 
that the element (1, 2) of B™ is given by 



I "- 1 

L 2 



x 



£ f l k W 2 + ±y) n - 1 - k i-y) k - 



k=0 



k 



The corresponding element of the RHS of Equation 1 is 



£ ? )(2 y r- k x k F k (x,y). 



fc=0 



k 



Then, using the Binet's form 

I n ~ 1 
L 2 



1 



sm^-Hs-^ 



£ ; (2y)»-*x*(a fc - i g 



Now 



\n— kk„ k 



m m ' - 



(2y + xaf 
a n (a - (3) r 



Analogously for the (3 part. After replacing n by 2n we get the identity 



rc-l 

k=0 



2n-l-k 



Now 



and 



BA 



(x 2 + 4y) n -*- 1 (-y)* = F 2n (x,y). (2) 



x 3 + 3xy x 2 + 2y 
x 2 y + 2y 2 xy 



tr(BA) = x 3 + 4xy = x(x 2 + Ay), 

|BA| = |Bj |A| = -y(x 2 y + Ay 2 ) = -y 2 (x 2 + %). 

Applying the previous method this time to BA we get that the element 
(1, 2) of (BA)" is given by 



(x 2 + 2v) ^ I ' \„n-l-2kfJ2 i / i„^n-l-2fe„,2fe/„2 

fc=0 V 



A: 



x 



(x 2 + %f- 1 - 2 V fc (^ 2 + 4y) A 



= (x 2 + 2y) £ 
fc=o 

On the other hand (see [1]) 



n — 1 — k 



„n-\-2k(Jl , / i„,\n-l-fc„,2fe 



x 2 + %)"^ i -V • 



BA = (2yl + xA)A = 2yA + xA 2 , 



so that 



(BA)' 



E ' (2y) n - fc A"- A; x fe A 2 



E ? (2y)»- fc x*A 



n—kk \ n+k 



fc=0 



fc 



Since the element (1, 2) of A n is F n (x, y), we have that the element (1, 2) 
of (BA) n is 



J2 I}) (2y) n - k x k F n+k (x, y) = -^- E \V\ (2y) n ' k x k (a n+k - (3 n+k ). 

k=0 V / fc=0 V / 



Now 



^ i (i) <«- w = ^ ± Q (*)-*(„)* 



;(2y + xa) n 
-a n (a-/?) n 



a — /3 
a" 



a — /3 
a 2n {a-j3) n 
a — (3 

In the same way 

' JT ( n X2yY~ k x k a n + k = J^-jr( n \2yy«- k (xP) k 



a-(3^ Q \kJ a - (3 ^ \k 

8 n 

Vv + xpy 



a — (3 



j " /3 n (-l) n (a-(3) r 



a — (3 
(3 2n {-l) n {a-(3) 
a — (3 

Subtracting we get 

/ a 2n [3 2n {-l) n \ 



(a-py 



a — (3 a — (3 



Replacing n with In and equating the two expressions for the element (1,2) 
we get 

(x 2 + 2y) £ ( 2U ~l " k \ x 2n - l ~ 2k (x 2 + ±y) 2n ^ k y 2k 
= (x 2 +4yTF 4n (x,y), 
that is 

(x 2 + 2y) £ f 2n "^ " ^ ^-^(x 2 + W 1 " V* = F 4n (x, y). (3) 



Using Equation 2 and replacing here n by 2n we obtain the identity 
(x 2 + 2y) J2 ( 2n ~l " k ) x 2n -^ 2 \x 2 + ^y) n ^ k y 2k 

= * 2 £ ( 4n ~l ~ k ) (* 2 + %) 2n " 1 - fc (-y) fe - 



3 Polynomials as Arguments 

First of all, from 

F n (x, y) = ?—%-, L n (x, y) = a n + /?", 

a — p 



solving for a and (3 we get 



Secondly 



It follows 



n = L n(x, y) + (a - (3)F n (x, y) 
an = L n{x, y) - (a - f3)F n (x, y) 



L 2 n (x,y) = a 2n + (3 2n + 2(-iry n 

= a 2n + /3 2n -2{-l) n+1 y n , 

(a-(3) 2 F 2 (x,y) = a 2n + (3 2n - 2(-l) n y n 

= a 2n + l3 2n + 2(-l) n+l y 2 . 

L 2 n (x, y) + (-l) n+ Hy n = (x 2 + 4y)F 2 (x, y). (6) 



Then 




a(L n (x,y), (-l) n+1 y n ) 


L n {x, y) + y/I%(x, y) + 4(-l)™+V 


2 




L n (x, y) + yjx 2 + 4yF n (x, y) 
2 
= a n {x, y). 



(7) 



p(L n (x,y), (-ir +1 y n )=p n (x,y). 

Then 

F n (L k (x, y), (-l) k+ Y) 

_ a n (L k {x, y), (-l) fc+ V) - P n (L k (x, y), (-W k+1 y k ) 

a (L k (x, y), (-1)*+V) - P (L k (x, y), (=1)*+V) 
__ a nfc (x, y) - (3 nk (x, y) 
a k {x, y) - (3 k (x, y) 
__ Fnkjx, y) 
Fk(x, y) ' 

As a consequence, we can see that F k (x, y) divides F kn (x, y), V/c, n. 
The analogous of Equation 7 for Lucas polynomials is 

L n (L k {x, y), (~l) k+1 y k ) = L nk (x, y). (8) 

Some examples. 

F 2n (x, y) = F 2 (x, y)F n (L 2 (x, y), -y 2 ) = xF n (x 2 + 2y, -y 2 ), 

F 3n (x, y) = F 3 (x, y)F n {L 3 {x, y), y 3 ) = {x 2 + y)F n {x 3 + 3xy, y 3 ), 
Fau(x, y) = F A {x, y)F n {L 4 {x, y), -y 4 ) = x(x 2 +2y)F n {x 4 +4x 2 y+2y 2 , -y 4 ). 

Equation 6 and Equation 8 can be used to derive many other identities. For 
instance the Simpson formula (see [1]) 

L n (x, y)L n+2 (x, y) - L 2 n+1 (x, y) = {-l) n y n {x 2 + Ay) (9) 

setting x = L k (x, y), y = {—l) k+1 y k becomes 

L n (L k {x, y), (-l) fe+ y) L n+2 (L k (x, y), {-l) k+l y k ) + 
-L 2 n+1 (L k (x, y), (-l) k+1 y k ) 

= (-ir({-l) k+1 y k T(L 2 (x, y) +A{-l) k+1 y k ) 
= (-y) nk (x 2 + Ay)F 2 (x,y), 



that is 

L kn (x, y)L k(n+2) (x, y) - L 2 Hn+1) (x, y) = (-y) nk (x 2 + Ay)F k 2 (x, y). 

Equation 2 with x = L k (x, y), y = (—l) k+l y k becomes 

2n - J - r\ , 2 ± _ r 2 (n-l-r), w ,rt 



L k (x,y)J2( (x 2 + %r 1 -^"- 1 ^(x, 2 /)(-y) 

r=0 V r / 

= F 2n (L fc (x, y), (-l) h+1 y k ). 
Now using the identity 

Ll(x, y) + 2(-ir+ 1 y n = L 2n (x, y), 
we have 

i^n (L fc (x, y), (-l) k+1 y k ) 

= L k (x, y)F n (L&x, y) + 2(-l) k+1 y k , _((-l)*+V) 2 ) 
= £fc(z, y)^n [L2k{x, y), -y 2k ) ■ 

Finally 



7 2n — 1 — r 



E r )(x 2 +iyr- 1 - r F 2 k ^ 1 ' r \x,y)(-yY k = F n (L 2k (x,y),-y 2k ) 

r=0 V r / 

Then, since — y 2k = ( — l) 2k+1 y 2k , Equation 7 allows to write 

E ( 2n ~ l ~ r ) (* 2 + ^T-^F^-^ix, y)(- y y k - F ^ X > y) 



r=0 



r / " ...... F2fc ^ y j 



Equation 7 and Equation 8 allow to derive new identities. For example, the 
identity 

yF n _i(x, y) + F n+1 (x, y) = L n (x, y) (10) 

upon substitution of x with L k (x, y) and y with (— l) fe+1 y fc becomes 

f u+i An-pQc, V) F k{n+1) (x, y) _ 

( } V F k (x,y) + F k (x,y) ~ ^*' Vh 

that is 

(-l) k+1 y k F k{n _ 1) (x, y) + F k{n+1) (x, y) = F k (x, y)L nk (x, y). 



Or the identity 

L l+2(x, y) + yL 2 n+1 {x, y) = (x 2 + 2y)L 2n+2 (x, y) + xyL 2n +i{x, y) 
that becomes 

£fc(n+ 2 )(^ y) + (-i) fc+ V£fc(n+i)(*> v) 

= L 2 k(x, y)L k{2 n+2){x, y) + (-l) k+1 y k L k (x, y)L k(2n+1) (x, y). 

Equation 6 can be rewritten as 

L 2 n (x, y) = (x 2 + 4y)F 2 (x, y) + (-l) n 4y n . 

Then 

a (y/x* + 4yF n (x, y), (-1)"/) = a n (x, y), 

(3 (7x 2 + %F n (x, y), (-l) n y n ) = -P n (x, y). 
It follows 

L k(2 n+l)(x, y) 

k(x, y) 
F2kn(x, y) 



/•2,, + J [ V/- -- WU-V. //)• (-l)V) = ^x y{ ' (11) 



:i2i 



^n (7x 2 + 4yF fc (x, y), (-l)V) = (a - /3) j _ . ? y) 
For the Lucas polynomials we have 

L2n+i (Vx 2 + 4yF fc (x, y), (-l)V) = (a - P)F k(2n+1) (x, y), (13) 

L 2n (y/x 2 + 4yF k (x, y), (-l)V) = W*> y). (14) 

Now identity 10 replacing x with ^/x 2 + 4yF k (x, y) and y with (—l) k y k 
becomes, using Equations 11, 12, 13, 14, 

{-l) k y k L k{2n „ l) {x, y) + L k{2n+1) (x, y) = L 2fcn (x, y)L k {x, y), 

and 

(-l) fc y fc F 2fc „(x, y) + F fc(2n+2) (x, y) = F fc(2n+1) (x, y)L k {x, y). 

In an analogous way the Simpson formula (Equation 9) becomes 

L 2 kn{x, y)L k(2n+2) {x, y)-{a- (3) 2 F k{2n+1) (x, y) = y 2nk L 2 k {x, y), 

and 

(a - /3) 2 i ? fc (2n-i)(^, y)Fk{2n+i){x, y) ~ L 2 2kn (x, y) = -y fc(2n ~ 1) L|(x, y). 
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